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Abstract 

We calculate the ground state energies of closed fc-strings in (2+l)-dimensional SU(N) 
gauge theories, for N = 4,5,6,8 and k = 2,3,4. From the dependence of the ground state 
energy on the string length, we infer that such fc-strings are described by an effective string 
theory that is in the same bosonic universality class (Nambu-Goto) as the fundamental string. 
When we compare the continuum fc-string tensions to the corresponding fundamental string 
tensions, we find that the ratios are close to, but typically 1% — 2% above, the Casimir 
scaling values favoured by some theoretical approaches. Fitting the A-dependence in a model- 
independent way favours an expansion in 1/N (as in Casimir scaling) rather than the 1/N 2 
that is suggested by naive colour counting. We also observe that the low-lying spectrum of k- 
string states falls into sectors that belong to particular irreducible representations of SU (N) , 
demonstrating that the dynamics of string binding knows about the full gauge group and not 
just about its centre. 



1 Introduction 



In this paper we calculate the ground state energy of a confining flux loop that winds once 
around a spatial torus of length I, with the flux in a higher representation than the funda- 
mental. We do so for a number of SU(iV) groups in D = 2 + 1 dimensions. This study 
complements recent calculations for fundamental flux loops [I] and the excited state spectra 
of these j2]. In a forthcoming paper we shall analyse the excited state spectrum of multiply 
wound flux loops [3] and similar calculations in D = 3 + 1 dimensions are under way. 

These numerical calculations have a number of motivations. Firstly there are old ideas that 
the relevant degrees of freedom of linearly confining SU(iV) gauge theories are string-like flux 
tubes and that these should be described by an effective string theory |3] which might describe 
all the important physics of the field theory, particularly as N — > oo |5J. Calculations such as 
ours can serve to test this programme; and, at the very least, to pin down the details of the 
effective string theory that describes the dynamics of long flux tubes [6J, [7] . Secondly, there has 
been dramatic recent progress in the construction of string duals to various supersymmetric 
field theories [8j. There is a great effort to extend this to gauge theories and to QCD (see 
[9] for recent reviews). Our calculations may provide useful information in the search for the 
appropriate construction. 

Our recent calculation [H [2] of the spectrum of a closed loop of fundamental flux as a 
function of its length I, demonstrated that this spectrum is accurately given by the simple 
Nambu-Goto string model [10] down to lengths much smaller than those at which an expansion 
in powers of 1/al 2 diverges. While this does not contradict the conventional approach of 
starting with a background consisting of a long 'straight' string and expanding in fluctuations 
around it [6] , it does imply that one can do considerably better by expanding instead in small 
corrections around the full Nambu-Goto solution. Hints of this possibility can already be 
discerned in some recent work, pj], [12], based on [6], [7] respectively, that showed that the 
next term beyond the Luscher term [6J in the expansion in powers of 1/al 2 is also universal 
and has the same coefficient as in the Nambu-Goto case. This fact, that we can identify the 
states as being essentially those of the Nambu-Goto model with the corresponding 'phonon' 
occupation numbers, means that the observed small deviations from Nambu-Goto have the 
potential to provide detailed and powerful constraints on the form of the additional string 
interactions that are needed. 

The /c-string that is the object of our study in this paper, can be thought of as a bound 
state of k fundamental strings. Earlier calculations have shown that for N < oo there do 
indeed exist such bound states both in D = 2 + 1 and in D = 3 + 1 [131 EH EE HE]. This 
binding cannot be readily incorporated into the usual analytic frameworks [51 [7J since the 
binding interaction presumably involves the exchange of small closed loops of string that are 
not under analytic control. Thus such strings are of particular interest. 

In the case of SU(iV) gauge theories in D = 2 + 1 there is an additional motivation. There 
has been significant recent progress towards an analytic solution of these theories using a 
variational Hamiltonian approach [17J. (See also related work in [IS]-) These calculations 
make predictions for both the fundamental string tension and for fc-string tensions. In [1] we 
found that the former prediction is very accurate and improves with increasing N, being only 
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about 1% off the correct value at N = oo. (Although this discrepancy is statistically very 
significant.) In this paper we shall provide a similar test for the predicted /c-string tensions. 

In the next section we describe our method for calculating the energy of a closed /c-string. 
We then present, in Section EJ our most accurate results for the /-dependence of this energy 
(which happens to be for the case k = 2) and compare to the expectation of the simple 
bosonic string model. In Section H] we calculate the /c-string tensions, extrapolate these to 
the continuum limit, and test the conjecture that they scale with the appropriate Casimir. 
We follow this, in Section [5j with a model-independent analysis of the iV-dependence of these 
string tensions. Finally, in Section El we look at the wave-functions of the ground and excited 
states and show that in fact they fall into representations of the group, rather than being 
determined simply by the center. A summary of this work has been presented in |19j . 



2 Methodology 

Consider a spatial torus of length / with all the other tori so large that they may be considered 
infinite. If the system is linearly confining, any flux around the torus will be confined to 'string- 
like' flux tubes. If the source from which such a flux would emanate transforms as ip — > z k ip 
under a gauge transformation, z, that belongs to the center, Zn, of the SU(iV) gauge group, 
the flux tube is called a /c-string. This is a useful categorisation since it is invariant under 
screening by gluons. For large / the ground state /c-string energy £7o,fc(0 gives us the /c-string 
tension, 

a k = hm — - — . (1) 

(In this paper o k is the tension of the ground state /c-string, <jf = o k= \ is the fundamental 
string tension, a-jz is the tension of a string carrying flux in representation 7Z, and a is a generic 
string tension. When there is no ambiguity, we drop the extra k label on Eq.) Whether stable 
/c-strings know about the full group or only about the center of the group, is one of the issues 
we address in this paper. 

As one varies / there will be a finite volume phase transition at l c = 1/T C , where T c is 
the value of the deconfining temperature. For / > l c the flux around the torus will be in a 
flux tube of length / and so we can ask what is the effective string theory that describes its 
spectrum. For / < l c there is no winding flux tube about which to ask such a question. 

The alternative approach of studying open strings using Wilson loops (or pairs of un- 
smeared Polyakov loops), is complicated by the contribution of the sources to the correspond- 
ing ground state energy. The closed string setup provides a theoretically cleaner framework 
for studying the effective string theory although for smaller N, where the phase transition is 
second order or weakly first order, the influence of the nearby critical point needs to be taken 
into account once I ~ l c . On the other hand open strings ending on point sources provide a 
way to address the important question of how perturbative physics at short distances evolves 
into non-perturbative physics at large distances, and whether there is a non-analyticity in 
this evolution as N — > oo. Comparing the SU(5) potential energy calculated in [20J with the 
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lower N potentials calculated in [21] points to the possibility of an interesting answer to this 
question. 

Possible effective string theories fall into universality classes. If we expand Eq(1) in powers 
of 1/al 2 around al, then for a bosonic string, where the only zero modes are those that arise 
from the D — 2 = 1 transverse translations, the first two correction terms are universal [HJ [12] 

More generally, if there are additional zero modes along the string, then we may have c ^ 1 
(although in that case only the first correction term has been shown to be universal). The 
simplest example of a bosonic string is the Nambu-Goto free string model, whose ground state 
energy is [10] 

which reproduces eqn(T5]) when expanded in powers of 1/al 2 . 

Our calculations are performed using standard lattice tachniques. We work on periodic 
hypercubic L x X Ly X L t lattices with lattice spacing a. The degrees of freedom are SU(iV) 
matrices, Ui, assigned to the links I of the lattice. The partition function is 

Z(l3) = ! Y[dU ie -^-^ u ^ ; P = — 2 (4) 

J j ag 

where U p is the ordered product of matrices around the boundary of the elementary square 
(plaquette) labelled by p and g 2 is the coupling, which in D = 2 + 1 has dimensions of mass. 
This is the standard Wilson plaquette action and the continuum limit is approached by tuning 
(3 = 2N/ag 2 — > oo. One expects that for large N physical masses will be proportional to the 
't Hooft coupling A = g 2 N, and this is indeed what one observes [22]. So if we vary (3 oc A^ 2 
then we keep the lattice spacing a fixed in physical units. 

We will consider a loop of flux that winds once around the x-torus, so that it is of length 
I = aL x . A generic operator <pi that couples to such a periodic flux loop is an ordered product 
of link matrices along a space-like curve that winds once around the x-torus, with the link 
matrices in the same representation as the flux. Correlations are taken in the t direction so 
that the energy of the loop is an eigenstate of the Hamiltonian (transfer matrix) defined on 
the xy space. Such a correlator may be expanded 

C(t) = (0j(t)0(O)> = Knimivac^e-^ ; Ei < E i+l (5) 

n=0 

where, if we have confinement, E Q is the energy of the lightest flux loop. Since the fluctuations 
that determine the error in the Monte Carlo calculation of C(t) may be expressed as a correla- 
tion function with a disconnected piece, the error is approximately independent of t while C(t) 
itself decreases exponentially with t. Thus one needs operators <pi that have a large projection 
onto the desired state, so that this state dominates C(t) at small t where the statistical noise 
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is still relatively small. This can be achieved using standard smearing/blocking/variational 
techniques (see e.g. [2"2|[T4"]). If we denote a Polyakov loop in the fundamental k = 1 represen- 
tation by l p , then typical operators for fc-strings will be <pi = Tr{lp~i}Tr{l p y and we include 
all such operators in our variational basis, including their smeared/blocked versions. 

Since E (l) becomes large for large I, it can only be accurately calculated for modest values 
of l\Jo and so we need to know what are the important finite / corrections to the linear term, 
E (l) ~ al if we are to extract a reliable value for a. This is even more important for heavier 
strings with k > 1, since these are heavier than the fundamental k = 1 string. In [1] it was 
shown that for the fundamental string, corrections to eqn([3]) are extremely small, even very 
close to l c and one can safely extrapolate using eqn([3]) from, say, /^/°7 ~ 3. However there is 
no guarantee that this will continue to be the case for k > 1. Indeed we expect that at fixed 
k the binding will vanish as iV —>■ oo, and so in that limit, at any fixed / however large, the 
lightest /c-string will actually be k fundamental strings, each satisfying eqn([2]) [23]. Moreover, 
the fact that all the strings must share the same critical length, / c , leads us to expect much 
larger deviations from eqn([3]) at small / for k > 1. In addition there have been suggestions 
[24] that k > 1 strings do not belong to the bosonic string universality class and have c > 1. 
Thus it is important to determine the /-dependence of the ground state A;-string and this is 
the issue to which we now turn. 

3 Central charge and / dependence of /c-strings 

We have performed finite volume studies for N = 2,3, 4, 5, 6, 8 at various values of the lattice 
spacing a. Here we shall present two of our most accurate studies for k = 2 strings. These 
are for SU(4) at (3 = 32.0 and SU(5) at (3 = 80.0. In units of the fundamental string tension, 
the lattice spacings correspond to a^/°7 — 0.2153 and a^af ~ 0.1299 respectively. 

In Table [D we present the k = 2 ground state energy as a function of / = aL x . We also 
present, for comparison, the corresponding energy of the fundamental k = 1 string. The 
critical value of / at which one loses the string is given by l c ~ 5a and l c ~ 8a for SU(4) 
and SU(5) respectively. Thus our range of closed string lengths extends down to nearly the 
minimimum possible value. 

There are many useful ways to analyse the /-dependence. Here we shall fit neighbouring 
values of / to a Nambu-Goto formula with an 'effective' central charge: 

m) = * (i - c *ff\jp) 2 (6) 

This will tell us what is the universality class of the effective string theory that describes k = 2 
flux tubes; in particular if lim^ooCe/j = 1 then this universality class is that of the simple 
bosonic string, just as it is for fundamental strings. In addition the rate of approach of c e //(/) 
to its asymptotic value will tell us something about the size of the corrections. 

In Fig. [T] we plot the value of c e ff versus / for the SU(4) calculation. We show values for 
the k = 2 string and also, for comparison, for the fundamental k = 1 string. We express the 
length / in units of the string tension o^. While this is the appropriate variable to use, as we 



4 



see from eqiiQ, it does mean that the scale in fixed physical units is slightly different for the 
k — 1 and k = 2 analyses. The horizontal error bars end at the two values of I used in the fits 
to eqn([6]). In Fig. [2] we have a corresponding plot for SU(5). 

We observe in Figs filTSl quite strong evidence for the claim that at long distances the k = 2 
flux tube behaves like a simple bosonic string, just like the fundamental flux tube. This 
contradicts the conjecture in [23] that lim^ooCg// = o^joj. (This latter ratio is, from our 
fits, about 1.35 and 1.52 for SU(4) and SU(5) respectively.) 

However it is clear that the corrections to Nambu-Goto are very much larger for k = 2 
than for k — 1. Indeed, if we limited our analysis to /i/cfc < 2.5 we might be led to agree 
with the conjecture of [23]. One way to quantify the size of the deviation is to add to eqn([3]) 
a correction term that is of higher order than the universal terms in eqn((2]), i.e. 

Fitting, for example, our SU(5) loop masses gives 



c' 



-0.09(2) : k = 1 
-0.93(5) : k = 2 



What is interesting about this result is that we see that it is not the k = 2 string that has 
an anomalously large non-universal correction, but rather it is the k = 1 string that has an 
anomalously small non- universal correction. A value d = 0(1) is precisely what one might 
expect for a generic effective string model in the Nambu-Goto universality class and that is 
what we find for the k = 2 flux tube. This observation appears to be both insensitive to the 
value of N and to the size of a: that is to say, it applies to the continuum limit. 

Although our calculations of Eq(1) for k > 2 are not accurate enough to provide evidence as 
compelling as for k = 2, they are consistent with k > 2 strings also belonging to the Nambu- 
Goto universality class, and so we shall assume that to be the case. Within our statisitical 
errors it also appears justified to assume that the dominant finite-/ corrections are those in 
eqn®. 



4 Continuum string tensions and Casimir scaling 

The prediction of [17] is that string tensions in D=2+l SU(iV) gauge theories should be 
proportional to the quadratic Casimir of the representation of the flux. This is, of course, the 
exact result in D=l+1, since in that case linear confinement is the Coulomb interaction. Based 
on the idea of effective dimensional reduction driven by a highly disordered vacuum, it was 
conjectured a long time ago [23] that this might also hold in D = 2 + 1 and D = 3 + 1. There 
is some additional evidence for this hypothesis from calculations of the potential between 
charges in various representations of SU(3) [26] . 
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For a given k the smallest Casimir arises for the totally antisymmetric representation, and 
this should therefore provide the ground state /c-string tension: 

= k{N - k) 

a f N-l ' 1 ' 

This is the part of the 'Casimir Scaling' hypothesis that we shall be mainly testing in this 
paper. For this purpose it is useful to have an alternative conjecture that possesses the correct 
general properties. A convenient and well-known example is provided by the trigonometric 
form 

(10) 



cr k _ sin N 



07 sin ? 

that was originally suggested on the basis of an M-theory approach to QCD [27]. 

In fact the full prediction of p2] for a k is more specific than eqnflHJ), since it also predicts 
a value for <jj in terms of g 2 , and including this gives: 

<r k k(N-k)(N+l) 1 



(g 2 N) 2 N 2 8tt 

Since we already know that the prediction for the k — 1 (fundamental) string tension is too 
high by about 1 — 2% (depending on ]V) p], testing eqn ffTTJ) and eqnflH]) does not come to 
exactly the same thing. 

We now turn to our calculations of o&. We will illustrate these using the example of SU(5) 
and then state our results for other values of N. 

The first step is to extract ground state energies from the Euclidean correlation functions 
that we obtain as a result of the variational calculation described in Section [2J The standard 
method is to fit the ground state correlator with a single exponential, |c| 2 exp{— E (l)t}, for 
t > to, and to find the minimum value of to that provides an acceptable fit. (Usually the upper 
limit of the range in t is chosen so as to exclude the very noisy values at large t.) The value 
of the ground state energy obtained by such single-exponential fits we label by S. Of course, 
since in general |c| 2 < 1, we know that there must be some contribution even at t > to from 
the excited states in eqnQ. In the range of t fitted this contribution will be masked by the 
statistical errors, but will lead to a systematic downward shift in all extracted values of Eq(1). 
Ideally this shift should be within the statistical errors but this is far from guaranteed given 
our inevitably imperfect control of the statistical analysis. To bound this systematic error, we 
also perform fits with two exponential terms where the energy of the excited state is chosen 
to be that of the first excited state (using S fits) with the same quantum numbers. This will 
usually provide an upper bound on the excited state correction to the ground state energy. 
We label the ground state energy obtained by such double-exponential fits by D. (Note that 
this choice will not necessarily bound the effect of such contributions to the ratio of string 
tensions.) 

We calculate the values of E (l) for strings of various k on lattices with spatial tori that 
satisfy > 3. We have seen in Section [3] that at such distances the Nambu-Goto expression 
in eqn([3]) accurately represents the /-dependence of the ground state string energy (within 
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statistical errors that are typical of the present calculation). We therefore use eqn([7j) to 
extract values of from the corresponding E (l). 

Having obtained values of at various values of (3 we now want to extrapolate the 
dimensionless ratios o^jof to the continuum limit. We do so using: 

— -j-L = —— + cia a s + c 2 a a f + ... 12 

where the Cj can be treated as constants over our range of /?. Typically extrapolations of this 
kind include just the 0(a 2 ) correction, and one drops the points at the largest values of a until 
a statistically acceptable fit is possible. We shall also perform a 0(a 4 ) fit so as to estimate 
the systematic error in the continuum extrapolation. 

We show the resulting values of cr k /af in SU(5) in FigsOEl where the calculations have 
been performed using single and double exponential fits respectively. We also show linear 
0(a 2 ) and quadratic 0(a 4 ) extrapolations to the continuum limit in each case. The resulting 
continuum values, together with the values of x 2 P er degree of freedom for the best fits, are 
listed in Table [2j Comparing Figj3] with FigJH we see that the S and D values are quite 
consistent although the errors of the latter are significantly larger. The corresponding contin- 
uum limits 9X6 . clS db result, also very similar. If we now compare the two different continuum 
extrapolations, we see that including an extra 0(a 4 ) correction reduces the continuum value 
by a little more than one standard deviation. In principle the fit with the extra correction 
term should be more reliable. However the coefficient of such a term is in practice largely 
determined by the calculations at the largest values of a, where there is always the danger 
that one is being influenced by the strong-to-weak coupling crossover where the power series 
expansion in a 2 breaks down. In any case, we see from this study that neither the use of 
double exponential fits, nor the inclusion of higher order terms in the continuum extrapola- 
tion, makes much difference to the final result. Our qualitative conclusion is that while the 
continuum ratio is close to the Casimir Scaling value of 1.5, there is a significant and robust 
discrepancy at the ~ 1.5% level. 

We have repeated the above calculation for SU(4), SU(6) and SU(8) and the resulting 
continuum values, using 0(a 2 ) extrapolations, are presented in Table [3J The first error is 
statistical and the second is intended to give some idea of the maximum possible systematic 
error, by looking at what one obtains from various combinations of S and D fits to <Jk and 
<jf. (As such it is not an estimate of the actual systematic error which would undoubtedly 
be significantly smaller.) Some remarks. Firstly, for SU(6) our best fit has a rather poor, 
but not impossible, x 2 ■ We have therefore doubled the errors so as to achieve a reasonable 
confidence level. This will ensure that the N = 6 value does not distort our analysis of the 
iV-dependence in the next Section. Second remark: for iV = 8 one of the two finite volume 
calculations has a very poor x 2 when we attempt to extract the string tensions. We therefore 
do not include the resulting values, that have very small errors but very small confidence 
levels, in our continuum extrapolations, but instead take from that study the string tension 
as calculated from the lattice whose size (in physical units) is the same as we use at other 
values of f3. (Some freedom as to which values to include at those (3 at which we perform finite 
volume studies, as well as some freedom in the fitting process, explains why the SU(5) values 
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in Table [3] do not coincide exactly with those in Table [2j) Note that as k increases, the loop 
energy also increases, and the errors become less reliable. Nonetheless, despite these caveats, 
we see from Table [3] that we can confidently conclude that while the string tension ratios are 
remarkably close to the Casimir Scaling values, and much further away from (MQCD) Sine 
Scaling, there is a definite discrepancy that is typically at the 1 — 2% level. 

It is interesting to note that if we take the actual values of ^/<jf/g 2 N as obtained in [1] 
together with the values of in eqn ljTTI) . we obtain the values of &k/&f listed in the last column 
of Table [3j We see that these values are largely consistent with our calculated numbers. That 
is to say, while the prediction of [17] for 07 is too high by ~ 1 — 2%, the prediction for ^Jokj g 2 N 
with k > 1 is actually consistent with the values we obtain. This is particularly significant for 
the k = 2 values which are the ones that we determine most accurately and reliably. 



5 TV-dependence 

An interesting feature of the Casimir Scaling hypothesis in eqnQ is that if we expand the 
expression for o^joj in powers of 1/N at fixed k, we see that the leading correction is 0(1 /N) 
rather than the 0(1/N 2 ) that one would expect from the usual colour counting rules. This 
observation [28] has generated some controversy [29j [16] . It is therefore interesting to see if we 
can learn something about the power of the leading correction from a more model-independent 
analysis of our results. 

We begin with our results for the k = 2 continuum string tensions, as these are more 
accurate than those for larger k. We attempt a fit with the ansatz 

&k—2 a b , 

k - z —2 n 31 

using the values p = 1 and p = 2, and show the best such fits in Figj5l 

For p = 1 our best fit gives a = 1.48(5) and b = 4.42(21), with x 2 /n d f ^ 2.7/2 = 1.35. 
Thus we infer that an 0(1/N) correction is quite consistent with our k = 2 calculations even 
if, as we see in Figj5l the full Casimir Scaling prediction is not. It is also interesting to note 
that this best fit is consistent with the trivial value of 0^=2/ / o r / = 1 for SU(3), even though this 
value has not been included in the fit. For p = 2, on the other hand, we find no acceptable 
fit at all, and in order to proceed we have to drop the N = 4 point from the fit. We now 
only have one degree of freedom, and although we obtain an acceptable fit, with a = 18.0(5) 
and b = —155(12), for this reduced dataset, we see from Figj5]that it has a wavy behaviour 
that is characteristic of an inappropriate functional form being forced onto data using large 
cancellations between different terms. This perception is reinforced when we note how much 
this fit deviates from the SU(4) value. We conclude from all this that our k = 2 results 
strongly disfavour a conventional 0(1 /N 2 ) correction. 

Complementary to the above large N limit, in which k is kept fixed, is the limit where k/N 
is held fixed. The extreme example is k = N/2; i.e. k = 2,3,4 for iV = 4, 6, 8 respectively. 
Of course the k = 4 string is quite massive, so it will have a large statistical error, and quite 
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possibly a large systematic error as well. Proceeding anyway, we fit with the ansatz 



with p — 1, 2. As we see in Fig® the best p = 1 fit is acceptable (with \ 2 — 1-5 for 1 degree 
of freedom) and gives a = 0.501(4) and b = 0.70(2). It is interesting to note that this value 
of a is consistent with what one obtains at N = oo from Casimir scaling, although, as we see 
in Figj6j the latter is actually far from the calculated ratios at finite N. The p = 2 best fit is 
noticeably worse, with x 2 / n df — 2.0 and giving a = 0.567(3) and b = 1.78(5), but it is not so 
bad that it can be completely excluded. Nonetheless it mildly reinforces our earlier conclusion 
that 0(1/ N 2 ) corrections are disfavoured. 

6 Excited strings : group or centre? 

As we have seen (here and in earlier work) there do indeed exist non-trivial fc-string bound 
states. The ground state fc-string is stable, with string tension au < kcjf at finite N, and we 
have provided strong evidence that it falls into the universality class of the simple Nambu-Goto 
string model. Thus one may conjecture that there will be a spectrum of Nambu-Goto-like 
excitations of this fc-string, at long distances, for which the scale is set by o^. 

We also found that the string tension of the ground state /c-string is remarkably close to 
the prediction of Casimir Scaling i.e. that it is proportional to the quadratic Casimir of the 
totally antisymmetric representation of k fundamentals. If this is more than mere coinci- 
dence, it implies that we should also expect to see some sign of flux tubes corresponding to 
other representations, with string tensions (approximately) proportional to the corresponding 
Casimirs. Of course, since these Casimirs are larger, such flux tubes will not be stable, but can 
be screened by gluons to the totally antisymmetric one. If the amplitude for such screening 
is large, such 'resonant' states will be 'broad', difficult to identify, and will presumably play 
no significant dynamical role. However if the amplitude for such screening is small enough, 
then it should be possible to identify such 'resonant' flux tubes in the excited state spectrum 
of the fc-strings. Indeed, in this case, one may ask if such a 'resonant' flux tube possesses its 
own spectrum of Nambu-Goto-like excitations. 

By contrast, in the free Nambu-Goto string model, which has been found to work very well 
for fundamental strings that wind once around a torus, [U [2] there are of course no bound 
state fc-strings; we simply have = kaf (as indeed one expects to have in the SU(iV) gauge 
theory at N = oo). One may speculate that the observed binding can be encoded through 
some correction terms to the Nambu-Goto action, leading to a bound state string for each 
value of k. Moreover it is possible that each of these stable bound states will be accompanied 
by its own tower of Nambu-Goto-like excited states, once the string is long enough. However, 
while there may be additional 'resonant' string excitations, it is not clear how these correction 
terms could embody the detailed group structure of the SU(iV) gauge theory in any natural 
fashion, and so one would not expect the excited states to be related in any way to the different 
representations of the SU(iV) gauge group. 




(14) 



9 



There is some evidence from calculations of potentials between sources in different repre- 
sentations that the corresponding fluxes do form flux tubes [26] , but such calculations do not 
go to large enough separations to be unambiguous. There is also some evidence from earlier 
calculations of closed loops, that the lowest excited k = 2 strings fall approximately into 
symmetric and anti-symmetric representations [13], although the observed near-degeneracies 
in these states [H] have been interpreted as saying that these are in fact the same states, 
and that the categorisation at a given k into different representions is illusory: that is to say, 
confining strings only know about the center of the group [3D] . 

To determine what is actually the case, one needs accurate calculations of a significant 
number of excited states at each k. This requires a large basis of operators enabling one to 
obtain good overlaps onto these excited states, paralleling the k = 1 calculations in [1]. Such 
a calculation is under way [3]. In the meantime, it is interesting to see what we can extract 
from the present calculation which has a rather limited basis, with overlaps on to the excited 
/c-string states that are only moderately good. 

We restrict ourselves here to the k = 2 excited state spectrum which is the one that we 
can determine most accurately. The operators whose correlators we have calculated are of the 
form Tr/ 2 and {Tr/} 2 where / is a Polyakov loop. There will be such operators corresponding 
to each of the smearing/blocking levels and this provides us with a non-trivial basis for our 
variational calculation of the k = 2 spectrum. This basis of operators allows us to construct 
Polyakov loops in the totally symmetric (2S) and antisymmetric (2A) representations: 



To construct k = 2 operators in other representations would require higher powers, e.g. 
Tr/ 3 Tr/^, which we have not calculated here. Thus we shall limit our analysis to the k = 2 A and 
k = 2S representations. Using the basic operators in eqn( |T5i) . we construct p = operators 
at time t by summing up Tt 2 a,2sI over the spatial coordinates. We call the resulting operators 
$2a(^) and $25 (t) where there is an additional suppressed label, b, for the blocking level of the 
gauge links used in the construction of /. These are our basis operators and our variational 
procedure will (ideally) produce the linear combinations that are closest to the actual energy 
eigenstates. 

One can immediately learn something interesting from these basis operators. Consider the 
normalised overlap of a symmetric operator <3>2S at blocking levels bs with an antisymmetric 
operator $ 2 yi at blocking levels b A : 



If these overlaps turn out to be large, then the screening/mixing effects are large, and there is 
not much point in trying to find different string states corresponding to different representa- 
tions. If they are small, then it becomes interesting to determine the representation content 
of the ground and excited /c-string states. 

In Table S] we list the values of these overlaps as calculated on a 24 2 32 lattice in SU(5), 
at a value of (3 that corresponds to a^crj ~ 0.13. This is both close to the continuum limit 



Tr 2A l = Tr/ 2 - {Tr/} 2 



Tr 2S Z = Tr/ 2 + {Tr/} 2 . 



(15) 




(16) 
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and ensures that the string is long enough, ly/ctf — 3.1. Our results show that the overlaps 
are remarkably small. Indeed they are all consistent with zero, within very small errors, 
except for those that involve the very largest blocking level. Blocked links at this level spread 
right across the spatial volume and we have convincing evidence from calculations on larger 
volumes that the somewhat enhanced overlap is primarily a finite volume effect. Nonetheless 
even these non-zero values are in fact very small. This provides striking evidence that the 
screening dynamics is a weak perturbation on the classification of states into representations 



We now turn to our variational estimates of the the k = 2 string eigenstates, that we 
obtain using the full basis of $2a(^) and &2s(t) operators. We project the resulting eigenstates 
onto the k = 2A and k = 2S subspaces. (Which we do by forming an orthogonal basis out of 
our non-orthogonal $2/1,25 operators.) We call the corresponding overlaps Oja and Ojs, where 
we expect |Oj^| 2 + |Ojs| 2 ~ 1 since the k = 2A and k = 2S subspaces are nearly orthogonal. 
If the subspaces were exactly orthogonal then it would of course immediately follow that the 
eigenstates would fall into exact k = 2A and k = 2S representations. What we find is that 
they appear to do so to quite a good approximation. An example is shown in FigJTJ Here we 
plot the overlaps (Oj^j 2 , |Ojs| 2 for the first five k = 2 states, J = 1, ...,5, as a function of the 
string length I, in a calculation in SU(6) at (3 = 59.4 (which corresponds to a^Jo f ~ 0.28). 
What we see is that the states are almost entirely either pure k = 2A or k = 2S except at 
certain values of I where the state changes from one representation to the other. Since this 
always happens to pairs of states at the same value of I and involves an opposite change in 
representation in this pair, it is clear what is happening: as we increase I the two energy levels 
cross at that I. So, for example, in Fig 13 the states that are ^ 2 and \I/ 3 for l\Jo ^ < 3 exchange 
their ordering in energy for ly/crf > 4. And a similar exchange occurs for ^4 and ^5, but at a 
larger value of I. What we see, therefore, is that the /c-string states belong to particular SU(iV) 
representations to a very good approximation. The apparent near-degenaracies observed in 
[T4"] are in fact accidental: they arise from the (natural) fact that the various excited states 
have energies that vary differently with I, so that they cross as / increases and when they do 
so we have degeneracies. As we have just seen, the first k = 2 crossing occurs at ly/o-f ~ 3.5 
which is precisely where such degeneracies have been previously noted to occur [T31 [T4] . It 
is unambiguous that what we are seeing here is not a single state, as suggested in [30J, but 
rather two that are nearly degenerate at appropriate values of I. 

If we focus on the lightest three states in FigJTl we see that the ground state is always totally 
anti-symmetric, while the first excited state is initially symmetric and then, for l\J(?f > 4, 
antisymmetric. This very different dependence on I between the ground state k = 2S state 
and the first excited k = 2A state is easily understood if we think of there being two separate 
Nambu-Goto towers of states labelled by the k = 2A and k = 2S representations respectively. 
In the Nambu-Goto model the energy levels would be given by 



where a-ji is the string tension of the flux in representation TZ and n is an integer that counts 
the 'phonon' excitations along the string. For the ground string state we have n = and 



of SU(iV). 




(17) 
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Eq(1) increases approximately linearly with / in the range of I relevant here. On the other 
hand the first excited state with n — 1 will clearly have a very different variation with /. 
This is illustrated in FigJS] where we plot the Nambu-Goto energy levels (divided by a/l) for 
towers of states in the k = 2 A and k = 2S representations of SU(6) using the string tensions 
that are predicted by Casimir Scaling. One sees how the very different /-dependence of the 
ground state k = 2S string and the first excited k = 2A string ensures that they cross for a 
value of l\/(Tf close to where we observe the crossing in FigJTl Our estimates of the actual 
energy levels are very roughly consistent with such a scenerio, although a reliable and precise 
comparison must await a dedicated study of the k > 1 energy spectrum [3]. In the meanwhile 
what we have is very good evidence that the fc-string states fall, to a good approximation, into 
representations of SU(iV), and strong indications, from the crossings of the states, that there 
exists a ground state flux tube in representations such as k = 2S where it is not completely 
stable. 



7 Conclusions 

In this paper we have studied /c-strings that wind around a spatial torus of length /. This 
setup is convenient because it involves no sources so one knows that, for / > 1/T C , all that one 
has is a closed string-like flux tube. It provides a clean way to determine the effective string 
theory describing the flux tube. Once this question has been addressed the second, equally 
important question of how the effective string theory at long distances matches onto short- 
distance asymptotic freedom, can be conveniently addressed through the calculation of the 
potential between sources as a function of the distance between these sources. For example, 
one can ask whether the running coupling in the 'potential scheme' 

V (f) = -C^M (18) 
r 

(C is the quadratic Casimir for the representation of the sources) develops a non-analyticity 
at some critical distance as iV — ► oo. 

The main conclusions of our study of such closed /c-strings are as follows. 

• We find that just like the fundamental k — 1 flux tube, k = 2 flux tubes are described 
by an effective string theory that is in the universality class of the Nambu-Goto model. If one 
parameterises the correction to Nambu-Goto by an appropriate power of 1/crl 2 , one finds that 
the coefficient is 0(1); that is to say, it is neither small nor large. This is in contrast to the 
case of the fundamental string, where such corrections are remarkably small pQ. 

• The k > 3 flux tubes are also consistent with being in the simple bosonic string univer- 
sality class. While the statistical errors are larger than for k = 2, they also exclude a central 
charge oc Uk/uj as suggested in [21]. 

• The ratio of the ground state string tensions, a^ja^ is within 1 — 2% of Casimir Scaling, 
but this discrepancy is statistically significant, and it is robust against systematic errors. On 
the other hand, the values of <7k/g 2 N are consistent with the predictions of [TTJ for k > 1; in 
particular for k = 2 where we have very precise results. From this point of view, the small 
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breakdown of Casimir Scaling that we observe can be entirely attributed to the fact that the 
predicted value of <7/ in [T7] is slightly higher than its actual value as obtained in [lj. One 
should perhaps not read too much significance into this agreement. The important point is 
that any discrepancy between the predictions for in [17] and the true values are no larger 
for k > 1 strings than they are for the fundamental string. 

• An analysis of how the k = 2 string tension varies with N strongly suggests that the 
corrections to the N = oo limit, <Jk=2 — 2oy, come in powers of 1/N (as in Casimir Scaling) 
rather than in powers of 1/N 2 , as suggested by standard colour counting arguments. There 
is also some statistically weaker evidence for this from our analysis of the iV-dependence of 
Ok=N/2 for N = 4,6, 8. 

• The fc-string eigenstates fall, very accurately, into representations of SU(iV). That is 
to say, fc-strings know about the full group representation of the flux, and not just about its 
transformation properties under the centre. The near-degeneracies previously observed be- 
tween some states assigned to different representations, should be interpreted as level crossing 
(as a function of the string length I). The intriguing possibility that stable fc-strings (and per- 
haps even unstable 'resonant' strings) may possess their own towers of Nambu-Goto states, 
receives some support from preliminary precision results on the fc-string spectrum [3]. 
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SU(4) = 


= 32.0 




SU(5) P = 


: 80.0 




aE (k = 1) 


aE (k = 2) 


L x 


aE (k = 1) 


aE (k = 2) 


6 


0.1637(8) 


0.2428(17) 


10 


0.1015(6) 


0.1720(20) 


8 


0.2976(6) 


0.4241(11) 


12 


0.1519(5) 


0.2485(12) 


10 


0.4075(5) 


0.5716(10) 


16 


0.2344(5) 


0.3747(9) 


12 


0.5105(7) 


0.7043(23) 


20 


0.3089(6) 


0.4867(11) 


16 


0.7073(13) 


0.9656(35) 


24 


0.3827(6) 


0.5938(15) 








32 


0.5242(11) 


0.8090(28) 



Table 1: The masses of the lightest k = 1 and k = 2 flux loops as a function of their length 
/ = aL x , for SU(4) and SU(5) at the indicated values of P = 2N/g 2 . 



\im a ^ Q a k=2 /(7f ; SU(5) 


fit 


S X 2 /n d f 


D X 2 /n d f 


0(a 2 ) 
0(a 4 ) 


1.528(4) 1.3 
1.520(5) 1.1 


1.524(7) 0.9 
1.509(11) 0.7 



Table 2: Continuum limits of the extrapolations in Fig j3l4l with quality of best fit. 



k 


N 


(Tk/cTf 


Casimir 


Sine 


Nair 


2 
2 
2 
2 


4 
5 
6 

8 


1.3552(22)© 
1.5275(26) ( 9 6 ) 
1.6234(66) ( 9 8 ) 
1.7524(51)0 


1.333.. 

1.5 

1.6 
1.714.. 


1.414.. 
1.618.. 
1.732.. 
1.848.. 


1.361 
1.529 
1.629 
1.741 


3 
3 


6 

8 


1.8522(48) (^ 4 ) 
2.174(19) (?) 


1.8 
2.143.. 


2.0 
2.414.. 


1.832 
2.177 


4 


8 


2.366(11) ft 1 ) 


2.286.. 


2.613.. 


2.322 



Table 3: fc-string tensions compared to Casimir and Sine Scaling conjectures in eqns (l9lfT0l) 
and the Nair prediction in eqn tflTl) . First error statistical, second estimate of systematics. 
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b A 


&5 


O45 


b A 


bs 


Oas 


b A 


bs 


Oas 


1 


1 


-0.0003(8) 


2 


1 


0.0006(6) 


3 


1 


0.0006(5) 


1 


2 


-0.0005(8) 


2 


2 


0.0005(7) 


3 


2 


0.0005(7) 


1 


3 


-0.0006(8) 


2 


3 


0.0002(8) 


3 


3 


-0.0001(10) 


1 


4 


0.0002(8) 


2 


4 


0.0006(10) 


3 


4 


0.0008(12) 


1 


5 


0.0025(6) 


2 


5 


0.0047(8) 


3 


5 


0.0066(11) 


4 


1 


0.0006(6) 


5 


1 


-0.0002(5) 








4 


2 


0.0005(9) 


5 


2 


0.0005(8) 








4 


3 


-0.0003(11) 


5 


3 


0.0009(9) 








4 


4 


0.0010(13) 


5 


4 


0.0038(12) 








4 


5 


0.0110(13) 


5 


5 


0.0290(15) 









Table 4: Overlaps of Polyakov loops in the k=2A and k=2S representations, with blocking 
levels bA and bs respectively. In SU(5) at f3 — 80.0 on a 24 2 32 lattice. 
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Ce ff u 
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Figure 1: The effective coefficient of the 7i/3al 2 term in the Nambu-Goto expression for the 
ground state energy, eqn([BJ, for the range of lengths indicated. For k — 1, •, and k — 2, o, 
strings in SU(4) at = 32.0. 
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Figure 2: As in FigfT]but for SU(5) at /3 = 80.0. 
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0.025 0.05 0.075 0.1 0.125 0.15 

a 2 (T f 



Figure 3: k = 2 string tension in SU(5), extracted using single exponential, S, fits, as a 
function of the lattice spacing. Continuum extrapolations linear and quadratic in a 2 are 
shown. 



a k=2 
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Figure 4: As in Figj3l but using double exponential, D, fits. 
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2 3 4 5 

h/cfor N= 6, k=2, aVo=0.27832 



Figure 7: Overlaps of the five lowest states in the k = 2 spectrum of SU(6) at /? = 59.40 
(a ~ 0.12 fm) vs. the string length. In blue(red) are the overlaps onto the antisymmet- 
ric(symmetric) representation. 



k=2,N=6 




Figure 8: The spectrum of energies vs. the string length from the Nambu-Goto model for 
SU(6) and k = 2. Blue(red) denotes the energies of the antisymmetric (symmetric) represen- 
tation. 
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